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^ '. Abstract 

We consider a band of fermions in two space dimensions with a flux phase 
' (relativistic) dispersion relation coupled to a local magnetic impurity via an 

s — d interaction. This model describes spinons of a flux phase and it is also 

: 

a qualitative model of the quasiparticles in a d x 2_ y 2 superconductor. We find 
' a zero-temperature phase transition at a finite coupling constant between a 

Ov 

' weak coupling unscreened impurity state and a strong coupling regime with 

a ; 

a Kondo effect. We use large- iV methods to study the phase transition in 

-d ■ 

this Kondo system away from marginality. The Kondo energy scales linearly 

° 

^ ■ with the distance to the transition . The zero-field magnetic suceptibility 

> : 

J^j ■ at zero temperature diverges linearly. Similar behavior is found in the T- 

matrix which shows a resonance at the Kondo scale. However, in addition to 
this simple scaling, we always find the presence of logarithmic corrections-to- 
scaling. Such behavior is typical of systems at an upper critical dimension. 
We derive an effective fermion model in one space dimension for this problem. 
Unlike the usual Kondo problem, this system has an intrinsic multichannel 
nature which follows from the spinor structure of 2 + f -dimensional relativistic 
fermions. 

PACS numbers: 71.27. +a,75.20.Hr 
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I. INTRODUCTION 



It is by now well understood that the presence of a small concentration of magnetic 
impurities into an otherwise non-magnetic metallic host can affect dramatically the low- 
temperature properties of the system. The prototype of these interesting phenomena is the 
Kondo Effect. PJ At the single- impurity level, there is a non-perturbative crossover between 
a Curie- Weiss law behavior (in which the impurity behaves like a non-interacting localized 
magnetic moment) at high temperature, and a strongly interacting regime in which the 
magnetic impurity and the band electrons form a singlet ground state. In other words, in 
this low-temperature (strong-coupling) regime, the band electrons conspire to screen out the 
spin (and magnetic moment) of the impurity. At zero temperature a similar crossover occurs 
as a function of an external magnetic field. This picture has been developed by the concerted 
use of the renormalization group "exact" numerical simulations and scaling ||, exact 

solutions via the Bethe-Ansatz |^,[| and large- iV expansions J?|||. 

From the point of view of scaling, the Kondo problem is a typical situation in which a 
trivial fixed point, which describes band electrons decoupled from the magnetic impurity, is 
destabilized by a marginally relevant perturbation, the coupling to the magnetic impurity. 
This leads to an asymptotically free renormalization group flow with a /5-function which 
is quadratic in the coupling constant. Marginal perturbations appear in critical systems 
at a critical dimension. The standard s — d Kondo Hamiltonian is effectively a model of 
one dimensional chiral fermions coupled to a single magnetic impurity through the forward 
scattering channel ||. Clearly, in this case we are at the lowest critical dimension. This 
is a direct consequence of the fact that the band electrons have a Fermi surface where the 
density of states is finite and essentially constant. Thus, the Kondo effect is ultimately due 
to the availability of states in the electron band which can efficiently screen the impurity 
spin no matter how weak the exchange coupling constant may be. 

Some time ago Withoff and Fradkin (WF) || considered a genaralization of the Kondo 
problem to systems in which the density of band electron states may actually go to zero at 
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the Fermi energy. They showed that if the density of states of the electron band vanishes 
at the Fermi energy as a positive power of the energy, the Kondo effect is suppressed for 
small values of the exchange constant and that the Kondo screening only happens beyond a 
critical value of this coupling. In fact, it is easy to see that the exponent r of the one particle 
density of states, N(E) oc \E — E F \ r plays a role here quite analogous to the distance to 
the lower critical dimension d — d c in critical phenomena. WF showed, using a combination 
of a "poor man's scaling argument" and a large- N limit, that at least for small values of 
the exponent r, this is the correct picture. Quite generally, if the exponent r > this is a 
non-marginal Kondo system. 

There are a number of systems of physical interest where this situation does arise. A sim- 



ple example are the fermionic excitations of a quantum antiferromagnet in a flux phase []T0 
More importantly, the normal state excitations of a <i-wave superconductor (with symmetry 
d x 2_ y 2) behave precisely in this fashion ||11|| . In the vicinity of each node of the gap function 
(hence, the use of the term "node" hereafter), the dispersion relation for the normal quasi- 
particles is linear in the momentum. Thus, sufficiently close to the node, the quasiparticles 
have an effective relativistic-like dispersion. In the theory of superconductivity (isotropic or 
not) [O the dynamics of the quasiparticles is usually pictured in terms of Nambu spinors. 



For the case of a d-wave superconductor, Nersesyan, Tsvelik and Wenger |L3[ have shown 
that this approach leads to effective Hamiltonian for the quasiparticles which takes the form 
of a massless Dirac Hamiltonian for each node of the gap, with the "speed of light" equal to 
the Fermi velocity. (Naturally, relativistic massless Dirac fermions themselves always have 
this property.) In this paper we consider a model which describes properly the coupling of 
flux-phase fermions to a magnetic impurity. It turns out that this model can also be used 
to describe the coupling of a local magnetic impurity to a d-wave superconductor, including 
pair-breaking effects. In a separate publication we will discuss in detail the problem of a 



magnetic impurity in a d-wave superconductor in more detail fll] , |T5| . 

In this paper we reconsider the Kondo effect for non-marginal systems. The model 
has fermions coupled locally to a magnetic impurity. The fermions are assumed to obey 



a relativistic-like dispersion law and hence a density of states vanishing linearly with the 
energy. For simplicity we consider models with just one species of "relativistic" fermions. 
We will refer to them as "having a single node". In particular we will discuss the case of 
an impurity coupled to an electron band with a density of states that vanishes linearly with 
the distance to Fermi energy Ep = 0. This case was not examined by WF who found that 
the singularity structure changed as soon as r > 1/2. We will show in this paper that at 
r = 1 the additional singularities conspire to give simple scaling laws modified by logarithmic 
corrections [|T6 1 . This picture is strongly reminiscent of a critical system at an upper critical 



dimension. 

As in the conventional Kondo problem, here too we can construct an effective one- 
dimensional theory. However, when one carefully reduces the 2 + 1-dimensional fermions 
with a relativistic-like energy-momentum dispersion to an effective 1 + 1-dimensional model, 
one finds that there are at least two angular momentum channels that are always coupled. 
Thus, these are all inherently multichannel Kondo systems. In particular, it is always the 
I = channel coupled either to £ = 1 channel or to I = — 1 channel. Which pair of angular 
momentum channels are actually coupled depends on the Parity of the node (or cone) to 
which the channels belong. We find two equivalent ways to represent the dimensionally 
reduced model. We have the freedom to choose the effective 1 + 1-dimensional fermions to 
have a local kinetic energy and, hence, to behave like conventional right movers. However, 
in this picture, the effective interaction with the impurity becomes non-local. The other, 
alternative picture, is to have a local coupling between right moving fermions which now 
have a non-local kinetic energy. In fact, due to phase space factors, the density of states of 
the effective 1 + 1-dimensional model goes to zero linearly with energy. The linearity of the 
density of states reflects the fact that in the original problem the fermions move in two-space 
dimensions. It is this feature what drives the system out of marginality and what generates 
a critical coupling constant below which the width of the Kondo resonance vanishes. 

In all cases of physical interest, Parity is an exact symmetry. This means that the 
number of nodes (cones) is even and that there should be as many cones with positive 



parity as there are with negative parity. In the case of both the <i-wave superconductors 
and of the flux-phase fermions, this property follows from the alternating signs of the gap 
function. In contrast, if either parity or time reversal were broken, all four cones would have 
the same properties under parity. However, in that case there would always be a non-zero 
energy gap on the entire Fermi surface. Thus, in general, the effective Kondo hamiltonians 
always have an exact degeneracy (angular momentum channels). The spin symmetry is 
the usual SU(2) spin rotation invariance (which here we call color). Here we will work 
with the SU(N C ) generalization of this symmetry, with the physical N c = 2. The angular 
momentum degeneracy leads to and SU (Nf) flavor degeneracy. For a problem with one 
node, Nf = 2. When more than one node is considered, Nf > 2. We will also consider 
impurity scattering amplitudes which may change the angular momentum of the fermions. 
The basic and simplest model is worked out in section |I[ 

In the second part of the paper we use large- N methods to investigate the behavior of 
these systems near the critical coupling |7|]. We find that there is still a Kondo scale once 
we go over the critical coupling. However, the position and the width of the resonance are 
not related anymore in the very simple way they are in the usual Kondo effect. The details 
are worked out in section [TTTj . The most salient feature of our results is the presence of 
logarithmic corrections to simple scaling in all quantities of physical interest, including the 
zero-field paramagnetic suceptibility (at T = 0), the Kondo scale and the T-matrix for bulk 
fermions. 

The large- N c theory that we present here suggests the following scenario. For small values 
of the coupling constant J, the impurity is free and effectively paramagnetic. In contrast, 
for large values of J, the impurity is screened by the fermions. The critical coupling J c is 
dimensionful and scales with the energy cutoff. However, the fact that these systems are 
inherently multichannel Kondo systems suggests that for J > J c the impurity is actually 
overscreened ]T7|]. In the language of the Renormalization Group this requires a flow with 
two finite fixed points: one infrared unstable fixed point for the paramagnetic-Kondo phase 
transition and an infrared stable to describe the multichannel behavior of the Kondo phase. 
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As in the conventional Kondo problem, the large- N c theory can only describe the formation 
of the Kondo singlet which in a non-marginal system can only take place at finite coupling. 
In marginal Kondo systems multichannel behavior is found as a next-to-leading order effect 
in the 1/N C expansion or to third order in perturbative scaling In the systems of 

physical interest mentioned above where this model applies, the number of channels Nf is 
always larger than the number of colors. In such a situation one can imagine that there 
exists a critical number of channels such that both fixed points actually coincide and below 



this critical number the transition dissapears altogether ||19|| . We should emphasize here 
that this physics cannot be accessed by a straightforward use of the large- N c expansion. 

The paper is organized as follows. In section [Q| we derive the effective 1 + 1-dimensional 
impurity models for systems of fermions in two space dimensions with relativistic-like dis- 



persion (flux phase). In section [TlJ we study these models in the large- N limit and determine 
their critical behavior. In section [TV| we calculate the propagator for the band fermions in 
the iV — > oo limit and use it to derive the T-matrix. In section [V] we consider a model with 
the most general type of scattering process for fermions with only one node coupled to a 
single impurity and determine the phase diagram. Section [VT] is devoted to the conclusions 
and to the discussions of the similarities and differences between the problem we discuss 
here and the conventional Kondo problem. In Appendix A we sketch the calculations of a 
few integrals. 



II. TOY MODEL: S-D HAMILTON! AN FOR A LINEAR DENSITY OF STATES 

In this section we want to study the s-d Hamiltonian a la Read and Newns 0, with the 
difference that we are not assuming a constant density of states for the conduction fermions. 
We are coupling an f- impurity to a bath of "relativistic" electrons, ie., electrons with a linear 
dispersion relation. We assume the chemical potential for the electrons (the Fermi energy) 
is zero. This is a very special case but it will turn out to be quite important for the case of 
a d-wave superconductor as it will be discussed in the following sections. 
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We start with a two-dimensional free fermion Hamiltonian with a linear spectrum and 
a Fermi velocity Vp which may correspond to the linearization of a fermion band of width 
2D, where D is an energy cutoff. We will consider a model of band fermions with N c = N 
spin components with N = 2 for the physically relevant case of spin £77(2). In general, this 
model will have a global SU(N C ) "color" (or spin) symmetry. This model has a well defined 
large- N limit. The free part of the Hamiltonian is 

Ho = J (^^(P) **(P) t 2 - 1 ) 

where ai and «2 are Pauli matrices, so that 

/ e~ id \ 

p- a = p 1 a 1 + p 2 a 2 = p \ (2.2) 

\e id / 

Eq.( |2.2|) can be diagonalized by expanding the fields into a linear combination of spinor 
eigenstates. Let's call them u±. This eigenmodes satisfy 

vfP ■ a u± = ±Vf\p\u± (2.3) 

A particular (and convenient) choice of eigenstates is given by 

ui ] = 4= and u { l ] = ±^= (2.4) 

V2 V2 

The fields \l/'s carry a spin index a apart from the spinor index % and j which we ignore at 
this point but we shall put it back in when dealing with the interaction with the impurity 
term. 

*i(p) = E uKp)Cx(p) (2.5) 

A=± 

so that 

^(p) = -L[ C+ (p) + c _(p)] 

^(p) = ^[C + (p)-Up)}e ld (2.6) 
Now we can expand the ( fields in a basis of angular momentum eigenmodes 



Cfc(p)= E ^C±,m(p) (2.7) 

m=— oo 

where m is an integer and ± indicates the positive/negative energy spinor eigenmodes of 
the Dirac equation given by Eq. (|2.1| ). We then have 

h = E r^v F \p\\(X, m {p)c+M-ciMC-M] (2-8) 

m=-oo J ° 2tt L 

In real space we have 

Url= /"^E E «^)^( m ,A(p) (2.9) 

Jo 2tt a=± Jo 2tt m ±^ 00 

To reduce to an equivalent one-dimensional problem we use polar coordinates where p = 
(p, 8), r = (r, 0) and p ■ f = pr cos(6* — 0). By making use of the expansion 



oo 

3 iz COS 1/3 



£ iV n (z)e in * (2.10) 



n=— oo 



where J n (<2) are the Bessel functions. They satisfy the relations 

i P rcos(0-0) + im g = imj ( w ym4> gince j_ = (_l)njj^ ( 2 .11) 

o 2n 



Using these in Eq.( |2.9j ) one gets 

^ l(r ,0) = f°^l g ^ m J m (pr)e-^[C m , + (p) + C m ,-(p)] 
Jo 2tt V2 ^.oo 

* 2 (r,^)= E ^ + V m+ ibr)e- 4{m+1) ^[C m , + (p)-C m ,-(p)] (2-12) 

Jo 2vr V2 

At the impurity site, i.e., for r — > 0, Jo(0) = 1 and J m (0) = for m ^ 0. So, we are left 
with 

*i(o)= r^4[co >+ (p)+co,-(p)] 

Jo Ztt a/2 

*a(0) = jf [U+(P) - C-i,-(p)] (2-13) 

This is telling us that the effect of having a cone-like dispersion relation for the fermions 
instead of a flat band will not only reduce the available density of states but it will also induce 
an angular momentum mode mixing. Even though we may regard the impurity as having a 
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(^-function spatial form factor as in the case of the usual Kondo effect, the coupling to the 
fermions will no longer be restricted only to the s-wave, or the spherically symmetric angular 
momentum channel. This feature will become evident in the following step. We want to 
rewrite the theory by introducing effective one-dimensional fermion operators. We use the 
radial component of incoming and outgoing waves. We find two "flavors" of right-movers 
given by 

ci(p) = vblC+,o(bl); c 2 (p) = Jjpf c+,-i(bl); for p>° 

(2-14) 

ci(p) = y\p\ C-o(bl); c 2 (p) = - yj\p\ C--i(bl); for p < o 

In other words, c\ are the fermion states with angular momentum I = while c 2 are the 
fermion states with angular momentum £ = — 1. We summarize these mappings as 

*i(o) ^/^fp/bUCp) (2-15) 

For a system with an SU(N C ) symmetry, the s-d interaction can be written in the form 

H imp = ^(0)^(0)^ = ^(0)7^^(0) f\ T* s fs (2.16) 

where flr^gfs is the impurity spin coupled to the electrons at the real-space position x — 0. 
For a system with an SU(N C ) symmetry, a = 1,...,N^ — 1. We will use the standard 
representation of spin operators (here, the generators of SU(N C ) defined in terms of the 
fermion operators / 7 as flr^fg. If the constraint f^f 1 = 1 is satisfied, the impurity is in the 
lowest (fundamental) representation of SU(N C ). Other representations can be constructed 
by changing the "filling fraction" of the impurity spin. Notice that, for SU(2), the only 
possible value of the filling fraction is equal to one. 

Once again we stress the fact that the angular momentum channels given by m = and 
m = — 1 couple to the impurity. This is an important difference with respect to the case of 
constant density of states, where only the m = channel couples to the impurity. We shall 
single out these two channels from the free part of the Hamiltonian to obtain an effective 
theory by replacing Eq. (|2.13|) into Eq. (|2.16|) 
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+ 



+ 9 - 

2 jo 



+ 9 - 

2 jo 



^~~~~~vp \ P \ [cUpXoAp) -Co t ,-(P)Co,-(p) 

o 2n 1 

°° v -^v F \ P \ [cii >+ (p)c-i,+(p)-cV(p)u-(p) 

Z7T 

p dp 



g dg 



2tt Jo 2tt 



d,+(p) Ta <b,+(g) + C!,-(p)r tt <b,-(g) 



+ Ci4p)T a (o,-(q) + Q,-(p)T a (oAq) 



s a 



p dp f°° q dq 



2tt 



2tt 



Cii,+ (p)r a C-i, + (g) + Cli,-(p)T*C-i,-(9) 



Cli, + (p)T a C-i,-(g) - Cli.-foKC-i.+fo) 



(2.17) 



By using the correspondence between the £'s and c's defined above we may write the Hamil- 
tonian in a more compact form. Since we found that this model contains effectively two 



"flavors" of right movers, we define a flavor index I 



, Nf with Nf = 2. The system 



has an effective SU(2) "flavor" symmetry which originates in the unavoidable mixing of 
angular momentum waves by the impurity. Similarly, we will write down the spin index 
a — 1, . . . , N c explicitly. The effective Hamiltonian is (repeated indices are summed) 



dp 
2^ 



\P\ 4Ap) 





r°° dq r— 








J-oo 2n " 


energ 


y is E(p) = 



9 1 C L(^) 



(2.18) 



shows that the coupling between the impurity spin and the band electrons is momentum 
dependent. 

We note here that there is an alternative representation of this model. It may be obtained 
by, instead of rescaling the fields by a factor of J]p\ as we have done here, defining a new 
momentum variable k = p 2 (an later extending k to (— oo, +oo)). In this form, the effective 
Hamiltonian becomes 



f°° dk^,.. t ... ... g f°° dk t .,, . f°° dk' + ... 

Heff = / —E(k)cl(k)cUk) + | / — ct(k) C / — cl(k' 



f a r a ap fp (2.19) 



In this case E(k) = sgn(/c) \/~k. In this representation, the interaction becomes local at the 
expense of a non-local kinetic energy. In fact these two representations are equivalent and 
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are the only ones compatible with local (anti)-commutation relations for the \l/-fields. This 
is the representation used by Withoff and Fradkin |J. In this paper we will not use this 
representation of the model. 

We may now proceed to study this theory with a functional integral formalism [^]. The 
Lagrangian, in imaginary time, is 



c(t)= I^ c I(p)(-^ + e (p)) C M 



+ 



Jo 

Nr. 



dp 
2^ 



\p\ 4*{p) 





r°° dq r— 


(TUJ 


J-oo 2tt v 



q\ c iM 



ft T a 
J a 'a/3 



+fl{^-h^U + e f (r) (flU-Q f ) 



(2.20) 



In Eq. (|2.20|) we have defined Jo = g(N c /2). We have also included a magnetic field h 
which we will choose to be a diagonal matrix h a and it will be defined below. The Lagrange 
multiplier field e/(r) has been introduced to enforce the constraint of charge (occupancy) Qf 
at the impurity site. In principle, the partition function separates into a sum of subsectors 
each of which is characterized by an impurity occupancy Qf. 

Following the standard large- iV c decoupling approach of Read and Newns [0, we write 
the spin operators in terms of fermions and find an effective four Fermi interaction. This 
interaction can be written in a simple form by making use of the well known identity, which 
holds for the generators of SU (N c ) 



N?-l 



2^ T L r a/3 - N c SapSua - 5 au 5 a p 
a=l 



(2.21) 



A Hubbard- Stratonovich (H-S) transformation is now introduced to decouple the fermionic 
quartic term which arises from this expansion. Up to an integration over the H-S fields v?;(t) 
and (fii{r), Eq. fl2.20|) is equivalent to 

dp 



£ 'M = /I^w(^ + £ W) C ^» + / 



oo 

-oo 27T 



|p| U{T)ftc la {p) + Vl {r)cl{p)U 



Nr. Nf 



■Jc 



d 



1=1 



fl [-^ - K + e f (r) \f a - 0/c/(t) 



(2.22) 



Eq. (|2.22|) can be rearranged, by field shifting and completing squares, in the form 
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C'(r) 



00 dp 
-00 2tt 



' d_ 
dr 



E{p) 







-1 



cUp)+[^ + E(p) \I\pMt) Mr) 



+ fl(r) 



d_ 



-1 



fa 



N, »' 



+ T £ El^(T)| 2 -g / e / (r) 



Jo 



(2.23) 



We now obtain an effective theory for the fields e/ and tpi by integrating out the fermions. 
The partition function is given by 

Z = J Dp Df DJ Dc De f Dip Dip* exp(- J dr£'(r)) 

= Z J D(pD(p*De f exp(-5 e// ) (2.24) 



where 



Zq = exp 



N c N f Tr I ^log(^ + E(p) 



is the partition function of free fermions. The impurity part of the effective action is 

N f 



s, 



eff 



- Tr log 

cr=l 



<9 ^ f dp 1 



w( r ) 



+ / rfr 



'N f 



Ji 



EN 



Vi=i 



Qf e f 



(2.25) 



Here we stress that the model of physical interest has N c = 2 and Nf = 2. 

The effective action of Eq. Q2.25|) has the standard form of reference @. The key difference 
here is the form of the free fermion Green's function which in this problem has a relativistic 
form. For the usual Kondo problem the magnetic impurity is coupled to system of band 
electrons with a constant density of states at the Fermi surface. In the model that we discuss 
here, the density of states of the effective fermions (the "right movers") is still constant 
but the interaction with the magnetic impurity has an explicit momentum dependence. 
This momentum dependence is such that the effective coupling at low momenta becomes 
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arbitrarily small. We will show below that up to a critical value of the coupling constant Jo 
there is no Kondo effect. This is in fact the result of reference 0. Notice that the momentum 
dependence of the interaction, is a direct consequence of the relativistic dispersion. Should 
any finite density of states arise, either by effects of a chemical potential or induced by 
disorder, a crossover to a conventional Kondo effect will occur. There is an important 
physical case in which a finite density of states is precluded by reasons of symmetry, and 
the Fermi energy has to be locked at zero. This is the case of the d-wave superconductors 



which we will discuss elsewhere 15 



We now consider the N c —>■ oo limit. Here, large N c means the limit in which the rank of 
the group of spin rotations becomes large instead of being 577(2). To proceed with the 1/N C 
expansion, we look first for static solution for y?/(r). From the SU(Nf) flavor symmetry of 
the effective Hamiltonian, there is a manifold of solutions which span the group SU(Nf). 
Clearly, all the solutions break SU(Nf) spontaneously. As in all impurity problems, it is 
impossible to break spontaneously a continuous symmetry of bulk fermions by coupling 
them to an impurity, which has a finite Hilbert space. Thus, this apparent spontaneous 
symmetry breaking is an artifact of the N c —* oo limit. In fact, we expect that it will 
already be restored by the leading 1/N C correction. This is precisely what happens in the 
large- N approach to the conventional Kondo problem [0. Thus, quantities which exhibit 
this apparent spontaneous symmetry breaking will get strongly corrected already in the next 
order in 1/N C . 

In what follows we will seek a static, symmetric, solution of the form <fi(r) = ip Q . We 
will now derive the form of the Saddle Point equations which will determine ip as a function 
of Jo and of the filling fraction of the impurity Qf. 

First we need to compute 



(2.26) 



+ p 

By working in imaginary frequency we automatically get the time-ordered expression form 
for Eq. ( |2.25D . To compute Eq. fl2.26j ) we introduce a lorentzian cutoff function f\(p) and 
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extend the integration over p to ±00. We have 

r°° dp \p\ ( A 2 

g Q {iu) = - / — 

J-ooZnp — iuj \p z + A 2 

We can go back to Eq. ( |2.25[) which now reads 
f°° du 



iuj 



A 2 



2ttv 2 A 2 - uj 2 



log 



U) 

A 2 



^2.27) 



S eff = -N c p — 1<* 

JO Z7T 



6j + ^(l-iV/|^o| 2 /V,A 2 ))' 



'-^#o| 
Jn 



-Qf e f 



[2.28) 



where 



1 



1 



lo § ( To 



2^1 - {uj/Kf \A 2 / 
and the magnetic field h has been set to zero. The saddle point equations are 

26/ 



(2.29) 



r°° da; 



+ Q/ 







(2.30) 



Jo 2tt e 2 + cj 2 (1 - ^l^ol 2 /(cj 2 ,A 2 )) 2 
where we should understand the integral as computed using some convenient adiabatic cutoff. 
The non-trivial solution for us>q is given by 



1 



™du 2u 2 f(ui 2 , A 2 ) (1 - N f \<p \ 2 f(uj 2 ,A 2 )) 
u 2tt 



5 2 + w 2 (1 _ Ar/ |^|2 f( u *,A*))' 



(2.31) 



It will prove useful to define the dimensionless variables x 
point equations Eq. (|2.30|) and Eq. Q2.31|) now read 

Qf , , f°° dx v 



and v 



The saddle 



TT = - s S n ( e /) / — 



Z/ 2 + X 2 



l _ A logx 
1— or 



^2.32) 



and 



1 

Jo 



A 



00 logx 
ax 



1-x 2 



1 - A 



log a; 
1 -x 2 



x 



+ * 2 (l-A^| 



(2.33) 



(ttvf) 2 jo 

where A = Nf\ip \ 2 /ttv 2 f . 

There exists a critical value for Jo, which we define as the value of Jo at the point where 
(fi departs from zero, for vanishing 6f It is given by 



1 



00 duo I uj 2 * 

n V ' I A 2 



A 



fix- 



log a; 



A 7T 2 



2.34) 



(ttvf) 2 Jo 1 — x 2 (tt-Uf) 2 4 

In the next section we use the SPE to extract the critical behavior of the system near J c . 
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III. SCALING AND ENERGY SCALES IN THE STATIC APPROXIMATION 



In the previous section we obtained the saddle-point equations (SPE) for this theory; 
however, these SPE have a singular behavior around the point v = 0, A = 0. In particular, 
an expansion in powers of A for small A is not possible for u — > 0. In this section we want 
to investigate in further detail the behavior of the SPE and the scaling behavior of A as we 
approach the critical point. We will find that, opposite to the situation in the usual Kondo 
problem, there are now two independent energy scales for A and v. 

We go back to the SPE expressed in their original form 



N c TT JO u 2 



Qf i 

sgn(e/) = - / dx ■ ^ ( 3A ) 

~ + x 2 (l - A M§ 



and 



1 f°° hgx ( logx \ x 2 

dx 1 12 1 - A i 12 77 ——72 (3-2) 



g Jo 1 - x 2 \ 1 ~ x J v 2 + x 2 (l - A ^ 



where A = Nf\ip \ 2 /nvp and g = AJ /(nv F ) 2 . 

In Eq. (|3.1| ), the integrand is the equal time propagator for the impurity, with the inte- 
gration variable x being the imaginary frequency scaled by the band electrons cutoff A. This 
equation sets up a scale for a crossover between two different behaviors, in much the same 
way as in the usual Kondo effect the imaginary part of the band electrons Green function 
sets the scale of the Kondo temperature. 

In the case of the usual Kondo effect, the expression equivalent to Eq.(|3.1|) is 



" % sgn(e/) = \ f dx ^TFTa? (3 ' 3) 

The case N c = 2 is, in a sense, a limiting case as can be seen from Eq. ( ^.3[ ). On one hand, 

the r. h. s. of this equation is a positive function. This fact forces sgn(ej) to be negative in 
order to have a solution. However, even in this case one gets 

Qf — = arctan ( - — -J — arctan (- — - J (3.4) 
2 \l e /l/ \l e /l/ 
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where A is some electron band cutoff. This suggests that for Qf = 0, the impurity level |e/| 
has to approach zero. However, if Qf = 1, there is no solution, since n/2 — arctan(A/|ej|) 
is a strictly positive number. Now, if N c is an integer bigger than 2, and letting A — ► oo, we 
find 

71 (\ ~ tf) = arctan ( A I ( 3 - 5 ) 



.2 Nj Vl e /i. 
Here we can distinguish two different regimes: Qf/N c — > 1/2 which corresponds to the 
case A << |e/|, and Qf/N c — > which corresponds to |e/| << A. 

However, for the system being discussed here, the situation is a little different and, 
actually, more complex since, as it turns out, now there are two different scales involved. In 
principle there is a scale set on x by A at the point where 

\ogx 1 logx 

l^-Ay— , or, (3.6) 

If the value of the frequency (and therefore of x) is small enough so that x 2 < 1, we can 
approximate 

x ~ exp ^— — ^ (3.7) 

For A small enough, the approximation is consistent. Working on this idea, we split Eq. ( [3.1| ) 

as 

Qf - 1 r 1/A dx « + I r dx 



.V K ./u V 2 + X 2 (l - A g^) 71 Je - 1/A V 2 + X 2 (l - A 



= A + / 2 (3.8) 

where we have used the fact that the only solution consistent with Qf/N c < 1/2 has e/ < 0. 

The SPE for A may be treated using a similar approach. Eq. ( |3.2| ) can be split into I[ 
and I2, where the first one is the corresponding integral up to e _1 / A , and the second one 
takes over from that point to infinity. Again we are interested in the small A and small v 
regime. The detailed computation of the integrals I±, I2, and I[, I' 2 is given in appendix [A]. 

We will be interested in the following limiting cases: 
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N c 7T 



In the regime in which v « e _1 / A << 1, with u, A << 1, the contribution from I2 
is neglegible and jf ~ ^ e x / A or v pa n e _1 / A << e _1 / A . Thus, in this regime, 
j^- « —. In this limit, the leading term from the other SPE is going to be (see 



Appendix [A], Eg. ( [ATI) ) 



7T 

T 



00, 



,-l/A 



7T 



1 

A 



(3.9) 



This gives 



A 



log 



1 

so 



log log 




■ 


4 90 


log 2 








it' 2 1_ 

4 90 





+ 



(3.10) 



and 



-l/A 

N c 



(si _ J_ 

U 90 



7T 



log 



(3.11) 



0/ 



We now consider the opposite regime, e « v « 1, where we obtain 



Of 

N 



1 

2 



7TU 



e- 3 / A and 



3 X i/ 3 

7T, 



,-l/A 



— » e- l ' A if ^ 
1/3 » Nc 



1 

2" 



(3.12) 



In this regime, clearly jf- — > \ and the other SPE gives 



1 

go 



7i 

T 



7T 



(3.13) 



We can get a solution for v by iteration on Eq.( |3.13j ), 
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Using Eq. ( |3.12| ) we get v which, after being replaced into Eq.( |3.13| ) gives a scaling 
form for A in this regime 



2 /V _ r 

I 4 90, 



7T fl Qf 



3 V2 N r 



1/3 



^lof 



7T 



2 ( - A 



(3.15) 



A 



Thus, 



log 



1 



»U 90 A 3 V2 AT C 



log log 


[1 (V _ ^ 

_7T V 4 90/ 


- -l 


log 


2 (ill _ 

7T V 4 90/ 


(i(i-3f) 


^ 1/3" 


-1 



(3.16) 

Notice that, for the physical case N c = 2, the two regimes Qf/N c << 1 and Qf ~ iV c /2 are 
identical. Only as N c grows large it becomes possible to distinguish one case from the other. 
However, in the case Qf ~ N c /2 the magnitude of A appears to go to zero as Qf — > N c /2~ 
even though the critical coupling is independent of the value of Qf. In contrast, in the case 
Qf/N c << 1 such unwanted feature is not present. It is worth to note here that, in spite 
of this apparent difficulty, we will show at the end of this section that the zero-temperature 
suceptibility at zero field has the same behavior in both regimes. Thus, this difficulty is not 
physically relevant. For simplicity, we will use the case Qf/N c « 1 to extrapolate to the 
physically meaningful case of iV c = 2. 

We can use these results to derive the /9-function for the coupling constant g$ in the limit 
N c — > oo. In order to do this we may use Eq. ( p.!3| ) for the regime in which jf- ~ ~, or 
Eq.( |3.11|) for the case jf — > 0. Starting from Eq. (|3.9| ) we may replace A in terms of v to 
obtain 

^ 2 1 



7T 



— — — ~ a v (log a 
4 go 



log v) 



(3.17) 



where a = njf- The /3-function keeps track of the flow of go as the cutoff A is decreased 
from very large values (infinite bandwidth). Hence, 
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(3(g ) = -A^ (3.18) 



In both limits Q f /N c -> and <2//iV c -> 1/2 we find 

P(go)= - go + -gl (3.19) 

9c 

where — = ^ . 

We immediately see that there are two fixed points: (a) a stable fixed point at go = 
and (b) an unstable fixed point at g = g c . The stable fixed point represents the weak 
coupling phase in which the fermions are decoupled from the impurity and the impurity 
spin is unscreened and there is no Kondo effect. The non-trivial fixed point at g = g c 
separates the weak coupling phase from a strong coupling phase that, in principle, should 
exhibit a Kondo effect. This unstable fixed point can be regarded as the usual marginally 
unstable fixed point of the standard Kondo problem, now pushed to a finite value of the 
coupling constant by effect of the reduction of the density of band states. Hence, in the 
present case, the interaction has to become strong enough so as to overcome the effect of the 
depletion in the number of states available in order to produce the Kondo screening. Thus, 
for systems with a small coupling to the impurity, there will be no Kondo screening and we 
should expect a decoupled impurity behavior. 

If go > g c , the system is in the Kondo screeening phase. We may define a Kondo scale 
Tk for the regime go > g c . From Eq.( |3.9| ) through Eq.( |3.1lD we see that, in terms of the 
physical parameters of the theory we have 



rr II N c 2 { Jo — Jc i , , 



>Q£\ 2 ( Jq 



(3.20) 



The results found here generalize the work of WF [[J. In that work the case of a density 
of states which vanishes linearly was not discussed on account that the structure of the 
singularities was not smoothly connected with the case of a constant density of states. The 
results of this section show that the main difference is the presence of logarithmic corrections 
to scaling in all the physical quantities. The case of a magnetic impurity coupled to a Fermi 
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system with a linear density of states seems to be analogous to the behavior of critical 
systems at their upper critical dimension. In contrast, the case of the constant density of 
states is the analog of a critical system at the lower critical dimension (marginal instability). 
The closest analog to this problem is the critical behavior of the Gross-Neveu Model in 3 + 1 
(space-time) dimensions ]2U| . 



The results that we find here hold for the simple regime of large N c which does not 
describe correctly the dynamics of the system with more than one fermion flavor. At this 
leading order in 1/N C , the existence of several channels is only reflected into a trivial degen- 
eracy and it does not lead to any new physics. However, as has been shown by Blandin and 
Nozieres |18[] , already the leading corrections in 1/N C will lead to non-trivial behavior in the 



"screened" phase, go > g c . In the standard Kondo problem, this behavior was confirmed by 
both the Bethe-Ansatz solution |17 and Conformal Field Theory ||21|| . We note here that, 



due to the non-local character of the effective one-dimensional theory discussed in section [H|, 
it is not possible to use neither the Bethe-Ansatz nor Conformal Field Theory to study the 
non-critical Kondo problem that we discuss here. Nevertheless, we expect that the physics 
of the phase with g > g c should be similar to that of the conventional overscreend Kondo 
problem. We also note that, in a more realistic model, one expects four nodes with two 
nodes of each chirality and, hence with more flavors and a more pronounced overscreened 
behavior. 

We close this section with a discussion of the zero-field susceptibility for this model. The 
coupling of a magnetic field only to the impurity spin has the form of a Zeeman term 

H h = -Y. f%U (3-21) 

<T=1 

Here, h a are the elements of a suitably chosen SU(N C ) diagonal generator, which also includes 
the appropriate normalization factors so as to set the Bohr magneton to one. We will need 
to ensure the tracelessness of the generator, and we can do so by keeping all the elements of 
order one with alternating signs. In this way we prevent the case of overpopulating one of 
the states that the Zeeman term splits with respect to the others, i.e., we explicitly avoid 
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cases like the one with N c — 1 elements being —1 and the last element being N c — 1. This is 
an important consideration to care about if we want to study the magnetic field crossover, 
i.e., if we were interested in obtaining the suceptibility at finite field. This was discussed in 



reference ||22|| . However, we will be interested only in the zero-field suceptibility and it will 
be sufficient to assume all the h a to be ±h, so that the sum of their squares gives a factor 
of N c . The magnetization M is given by 

Nc ~ Nc r°° rliii h 

a=X a=l J -°° 

The susceptibility at zero field is 



, f -iu{l-N f \ytff{u\A?))-K 



(3.22) 



X\ 



dM 



h=0 



dh 



h=0 



Nc (2 
A | v 



d 



dx 



3v J JO 71 v 2 



X" 



(3.23) 



Once again, we consider two regimes: Q//N c — > and Qf/N c — > 1/2. We discuss only the 
regime go > g c since the suceptibility is infinite below g c , where the impurity is unscreened. 



a) 



9l 

N c 



<C 1: 



For Of ss e 1 ^, we get 



T K KVp N c 



b 



9l ~ 1. 

N c ~ T 

For %f m | e~ 3/A , we have 



JcJo 
J() — Jr 



log 



(2Q/Y 



\ N C J \Jo-J, 



Jo 



(3.24) 



x(o) 



JcJo 



2nvp \Jo — J c 



log 



Jo 



7T V Jo — Jr 



1-511-91 

' Nr 



(3.25) 



Once again, we find logarithmic corrections to scaling. 

We see from Eq.( [3.25|) and Eq. (|3.24j ) that the zero- field susceptibility x(0), at zero tem- 
perature, has a behavior consistent with the general picture described above. It is finite 
for g > g c , but diverges as g approaches g c . This is consistent with the impurity being 
screeened (Kondo effect) for a value of the coupling constant bigger than the critical value. 
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However, as we approach the critical coupling constant from the screening (Kondo) phase, 
the susceptibility must diverge, since for values of the coupling lower than the critical value, 
A vanishes and the impurity effectively decouples from the band. In this regime one ex- 
pects the impurity will behave like a free spin in a magnetic field and, since we are at zero 
temperature, the susceptibility should diverge. 



In this section we derive the propagator for the band electrons to order N — > oo. We will 
use this expression to calculate the T-matrix in N — > oo limit. We will consider processes 
in which an electron with a given energy is initially in a state with well defined angular 
momentum, and it is scattered by the impurity into a state with possibly a different angular 
momentum. Thus, we will parametrize the T-matrix by the magnitude of the incoming and 
the outgoing momenta (namely, the energy of the state) and by the angular momenta of 
channels involved in the scattering process. This description is natural since at most only a 
pair of angular moementum states (£ = 0,-1 with our choice of node) are actually mixed 
by the impurity. 



In contrast, the computation of the correlation function (tyj(f) ^(r ')V f° r ^* an d r ' 
away from the impurity will involve a mixing of all the angular momentum channels. How- 
ever, in principle, we could add a source term and also decompose the source field in angular 
momentum modes, and keep only the modes that get mixed by the interaction with the 
impurity. In other words, we are considering the scattering of band electrons among angular 
momentum channels produced by the interaction with the impurity. In the case of a cou- 
pling to a flat band (usual Kondo problem), the only angular momentum channel that gets 
involved is the s — wave. Thus, in practice, we only need to consider the partial waves that 
actually get mixed by scattereing from the impurity. 

Therefore, it will be sufficient to work with the equivalent 1-dimensional fermions c\ and 
C2 while keeping in mind the correspondences defined in section [Dl The relevant part of the 



IV. T-MATRIX AND ONE PARTICLE GREEN FUNCTION 
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action is 



r°° duo f°° dp t ... 
Sf= / 7T" / 7T C ^ U ) + E(p))c la (p,uo) 

J-oo Z7T J-oo Z7T 



+ 



-oo Z7T J-oo Z7T J-oo ZTT k 

+ iptiu - u}')c] a (p,u)f a (u')} 



(4.1) 



We add sources r]j a (p, uo) and r)i a (p, oo) for the fields ci a (p, uo) and c\ a by adding the following 
source term to the action 

duo f°° dp 



Loo 27T Loo 2lt U ) + 4(P> ^)%r(p\ ^)) 



(4.2) 



After completing squares and integrating over the band fermion fields, we are left with 
(appart from a normalization factor) an action which depends only on the impurity fermion 
fields and the sources, 
N ° [duo 

cr=l ' 



. r 00 do; r 00 dp 

a u) J-oo 27T J-oo 27T 



duo f 00 dp + , . I 1 

-oo 27T J-oo 27T y — + E[J>) 



(4.3) 



After a few fairly straightforward manipulations and another integration now over the im- 
purity fermion fields, we can re-write the action as 

Nc -duo r duo' 



(4.4) 



where we have introduced the operators J CT 's and the kernel K which are given by 

\p\ 



f duo' f dp v 



i=i ■ 



and 



K{uo,uo') = (-ioo + e f ) 5(uo - oo') - J g J 



+ E(p) 
dp f doo" \p\ 



(4.5) 



N f 



2tt -iuo' + E(p) fr{ 



(4.6) 
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The generating functional for the Green's functions (at zero temperature) can be written as 

= exp(- Sffarf)) (4.7) 



so that correlation functions of the band fermion fields can be computed by differentiation 
over Z F 

(dfe «) C«, O)) = ( ^ Z F [,, ,t] )^ (4.8, 

We restrict the calcultation to the static solution case where ipi(fl — u) = (pi 5{VL — uj). Then 

where the T-matrix Tu>(Q) is given by 

MQ) = ^ jt (4-10) 

-iQ + e f + go(tQ) Ej=i N 2 



where go(^) was given in Eq. ( |2.26| ) and Eq. ( |2.27| ). 



The T-matrix of eq. ( |4.10|) exhibits several important features: 

a ) Eq. (|4.10|) gives the T-matrix at imaginary frequency. The T-matrix for real frequency 
(at zero temperature) is found by an analytic continuation to the real axis, iQ — ► Q. 

b ) If (pi and (p*r are non-zero, there exists a finite matrix element for the scattering of 
a band electron from a state with angular momentum I into a state with angular 
momentum I' (and viceversa). In contrast, for a flat band, the T-matrix is diagonal 
in angular momentum states. This matrix element is not invariant under SU(Nf) 
since it mixes the angular momentum channels. As we pointed out above, this is a 
consequance of the unphysical spontaneous breaking of the flavor (channel) symmetry. 
Thus only the trace of this T-matrix is physical. All other matrix elements will be 
suppressed by infrared divergent corrections to next order in the 1/N C expansion. 
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From the analytic continuation of the trace of the T-matrix one can find the position 
of the Kondo resonance Qk, the width Wk and the phase shift 5(Q). Using the fact 
that Yluli Wi\ 2 = Nf\(po\ 2 , we can write the trace in the form 

1TV 2 F A 



trT(fl) = T u (n) 



Q+ \e f \ + QA\og± -imA 



(4.11) 



which is valid for frequencies low compared to the bandwidth O C A but comparable 
to the dynamically generated scales. The phase shift 5(Q) is 



6(0.) = tan" 1 



0+ |e/| + QAlog 



(4.12) 



The imaginary part of the T-matrix is 



Im trT(O) 



\{nAv F ) 2 VL 



(4.13) 



(| e/ |+0(l-Al 0g M)) +(fAn) 

The function given by Eq. ( |4.13 ) has a peak at the Kondo resonance Qk, which is the 
solution of the equation 



e/ = n(l-AlogH 



One easily obtains the asymptotic solution 

N, /J - JA I 



(4.14) 



n 



K 



-7TV F - 



Qf V JqJ c 



log 



2 / ^Qj_\ 3 

IT \ N c J \J -Jc 



Jo 



+ 



(4.15) 



where we have kept only the leading logarithmic corrections. We see immediately 
that we have a resonance centered at fix, which is basically the Kondo scale ej that 



we found in section [II | but modified by the ubiquitous logarithmic corrections. In 
particular, the position of the Kondo resonance scales with the distance to the critical 
coupling constant and it goes to zero (the Fermi energy) at J c . 



The resonance width Wk can be read off from Eq. Q4.13 ) 



, 7rA 

w K *\n K \ — 



(4.16) 
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From Eq. ( |4.12j ) it is clear that the dimensionless parameter A = Nf\ipo\ 2 /wvp sets the 
scale for the range of the scattering amplitude. Thus, as far as transport properties 
such as the resistivity are concerned, A is the important parameter, whereas thermo- 
dynamic properties such as the suceptibility, depend entirely on the impurity energy 
scale ej, i.e., the Kondo scale Tr. The analysis done in the last section shows that, 
due to the strong logarithmic corrections, these two scales are completely different. 
Once again, in the conventional Kondo effect, Tk controls both effects. 

Hence, in contrast with conventional Kondo behavior, the interaction with the magnetic 
impurity induces a non-trivial renormalization of the propagator of the band fermions. This 
will happen even for J < J c , although this effect only occurs in higher order in 1/N C . Finally 
we stress that, due to the SU(Nf) flavor symmetry, there cannot be any channel mixing 
scattering processes in this model. However, in the following section, we will present a model 
which describes the most general impurity scattering processes for systems with only one 
node. In that model the channel (or flavor) symmetry SU (Nf) is broken explicitly already 
at the level of the impurity Hamiltonian. Thus we expect to find processes which will mix 
the various channels. 

V. GENERALIZED IMPURITY 

In the past sections we considered a situation in which the fermions coupled to the 
impurity spin only through their own spin density ^r^. Scattering processes of this type 
have the simplifying feature that particles and holes interact with the impurity independently 
from each other and exactly in the same way. However, in practice, more general scattering 
processes will be present. We also made the assumption that only one node of the two- 
dimensional fermions is present. 

Here, we will consider all possible scattering processes involving only one node. The 
coupling through the spin density has the very special feature that it is diagonal in the 
components of the spinors, i.e., the impurity does not mix particles with holes. In terms 
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of the effective Hamiltonian of Eq. ( |2.19 ), this diagonal coupling implies that there is no 



explicit mixing of angular momentum channels. Hence, the model has an SU(Nf) symmetry 
(with Nf = 2). Any process which mixes particles and holes will break the flavor symmetry 
explicitly. We will include now these processes. 

We will also incorporate the correct node structure. For example, in the case of a 
flux phase say, on a square lattice, there will be four different nodes corresponding to the 
symmetry points (±7r/2ao, ±7r/2ao) of the first Brillouin zone. As we discussed in section 
|, a (i-wave superconductor has an analogous node structure, at the points where the gap 
has nodes. We will discuss of this very interesting case in a separate publication fl5 |. It is 



straightforward to see that the main effect of a more general node structure is to increase 
the number of flavors in the effective model. In the one-node model of the previous section 
we found that there were two flavors and an SU(2) flavor symmetry. The existence of 
two flavors can be traced back to the spinor structure of the original problem. When the 
two-dimensional fermions have more than one node, the number of flavors becomes twice the 
number of nodes, Nf = 2N no d es - Thus, for a flux phase we will have Nf = 8 flavors. However, 
flux phases are even under parity. Hence, there are two pairs of nodes with opposite parity. 
The parity of the node is given by the relative sign of the two terms in the free fermion 
Hamiltonian of Eq. (|2.1|) . Thus, a change of parity is equivalent to the mapping 9 — > —6, 
in other words to a reflection across the X\ axis. In Section |J we showed that, with the 
choice of parity we made there, the channels with angular momenta £ = and £ = — 1 
get mixed by the impurity. Hence, if nodes with the other parity are also included, we 
will also get mixing between £ = and £ = +1, but no direct mixing between £ = +1 
and £ = —1. In Section || we also found that only the mapping between the components 
of the original fermion and the fermions of the effective one-dimensional theory carried 
information about which angular momentum components are mixed. This feature is always 
present. Thus, the net effect of having a multinode structure is to increase the number of 
flavors. Consequently, although the thermodynamic properties will be insensitive to parity 
assignments, the scattering amplitudes of the original fermions with the impurity will carry 
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the information of the parity of the nodes. This feature can be used to determine the nature 
of a superconducting state. 

In what follows, we will consider a system with four nodes, two with positive parity and 
two with negative parity, interacting with a single magnetic impurity. For simplicity we 
will assume that the impurity form factor is strongly peaked at small momentum, so that 
inter-node scattering processes can be ignored. We will keep, however, processes in which 
particles and holes may interact differently with the impurity and/or mix with each other. 
We parametrize these processes with four coupling constants Jj, with i = 0,1,2,3. The 
generalized impurity Hamiltonian is given by 

H imp = £J, *L(0) Ti r ^,(0) S (5.1) 

l,i 

where the indices I labels the nodes and a and (3 the spin components. To = I is the 2x2 
identity matrix and T, = <7j (i = 1, 2, 3) are the three Pauli matrices. Here Jo is the coupling 
constant that was used in the model of Section [H]. 
The free fermion Hamiltonian for each node 
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Ho = J ^{pM-^) (5.2) 

is invariant under the symmetry transformation C^f{p) — > a^(— p). The effect of this 
transformation on on H imp is to reverse the sign of both J\ and J%. This implies that the 
sign of J\ and J2 is irrelevant to the properties of the theory, which should depend only on 
the absolute values of these two coupling constants. In the expressions for the Fermi field 

*«(0) = ^/_"^ >/Wc«(p) and S a = fftjs (5.3) 

we have I = 1, . . . , N no d es and i — 1, 2. In what follows we use a single flavor index i = 
1, . . . , Nf with Nf = 2N noc i es . However, we will keep in mind that the impurity does not mix 
different nodes but it does mix the flavor components associated with the same node. 

The model of equivalent one-dimensional right-moving excitations has an SU(2) spin 
(color) symmetry and an SU(2N no d es ) flavor symmetry (associated with the channels) broken 
down to SU(N no de S )- The free fermion kinetic term has the form 
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Ho = J pv F ^c\ a (p)ci a (p) (5.4) 

ia 

and an impurity interaction term with the following structure 

2^ J 2^vlPll9l/^(?) c i9<(p)/j9 T « ( 5 - 5 ) 

where Ty is the coupling matrix which has a block diagonal form. For i,j associated with 
the same node T^- has the form 

T l3 = iV c (j <% + + J 2 ai j + J 3 a%) (5.6) 

In other terms the matrix has the form T ® I where T is the 2x2 matrix of Eq. ( |5.6| ) and 
I is the N nodes x N nodes identity matrix. We have also let the spin indices a, (3 to run from 
1, . . . , N c and used the identity of Eq. (|2.21|) . The form of the impurity Hamiltonian Eq. fl5.5|) 
shows that, if the coupling constants Jj are all different, the SU(2) flavor symmetry of each 
node is broken by the interactions but the symmetry involving different nodes remains intact. 

We will now proceed as in the previous sections and solve this model in the limit of 
N c — > oo. After a Hubbard-Stratonovich transformation, the quartic term of the euclidean 
action becomes (repeated indices are summed) 

dp 



J drati^M^ir) + J dr jT ± J\p] (a t {r) ft(r)c ai {p, r) + a*(r) clip, r)f a (r)) (5.7) 

where M a p = (Tq^) -1 . Since T is of the form T® I, then M has the same form, i.e., M ®I. 
In the following, M stands for the 2x2 matrix of Eq. ( |5.7| ). It can be easily shown that 



M 



N, 



c 



( ( ■ \\ 



9o-9-9 



N c M (5.8) 



y - (0i + m) go + 03 j 
where J» = gi/N% and g = (gi, g 2 , 03). In the case considered in section [H], i.e., when Jj = 0, 
Jo 7^ 0, we have a full U(2) flavor symmetry between the a-fields. In other words, the two 
channels play exactly the same role. When the coupling constants are all different, this 
symmetry is broken down to a U(l) x U(l), where one of the U(l) symmetries is generated 
by the identity and the other by M. 
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The N c — > oo limit is taken in the standard fashion. The Hubbard-Stratonovich fields 
Oi will be chosen to be an arbitrary vector for each node and the same vector for all nodes. 
Let V be a unitary transformation which diagonalizes the 2x2 matrix M, and let a = Va, 
so that | cr x | 2 -|- | cr 2 1 2 = 1 5"i | 2 + 1 5^2 1 2 - The eigenvalues of M can be rewritten as m±5m, where 

9o j s y/W 



m 



90-9 2 



and 



Sm 



The modified effective action can be written as 



-Nr. I — lof 
3 lH 



gl-g 2 



(5.9) 



-iu + e f + -N nodes ( |o"i| + |cr 2 | ) G (iu / 



+ -^N no d es N c 

r°° du ^ 



5m 



(|o"l| 2 + |cr 2 

g + e/ + -Nnodes (l^il 2 + |<5- 2 1 2 ) G (^)^ 



l^ 2 | 2 



+ ~N nod esN c [ ^ ■= + '7'/^ ) - <2/ e / 

2 \9o-V9 2 9o + V9 J 



The new Saddle Point Equations are 

1,00 du 



and 



Qj_ 

N r . 



0\ 



oo 2n -iu + e f + \N nodes (\ai\ 2 + |a 2 | 2 ) G (iuj) 



and 



(To 



00 dw 
01 I 7T- 



Gn(iu) 



-oo 27r -ia> + e f + \N nodes (|<ti| 2 + |5- 2 | 2 ) G (iu) 



go + VlF ° 2 J-oo2-k -iu + e f + \N nodes {\<?i\ 2 + \a 2 \ 2 ) G (iuj) 
We now solve these new SPE's. The solution will be cast in the form of a phase diagram 
which can be plotted in a y/cp — go plane (see figure [l]). In principle, we find three different 
solutions of the SPE's: 

a ) a i = a 2 — 0. This is the region of the phase diagram below the line go + \fW — 9 c- 
This is the weak-coupling phase. The magnetic impurity is effectively decoupled from 
the band electrons. To leading order in 1/N C the impurity does not interact with the 
band fermions and behaves like a free magnetic moment. Consequently the impurity 
spin suceptibility is infinite, there is no resonance and there is no Kondo screening. 



00 duo 



Go(iuj) 



(5.10) 



(5.11) 



(5.12) 



(5.13) 
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b ) Both a i ^ and a 2 7^ 0. There is no consistent solution of this form unless g 2 = 0. 
In other words, this case is possible only if the only non-vanishing coupling constant 
is Jq, with Jo > J C i in which case we do have the U(2) symmetry. This is the line on 
the axis go of the figure [1| for g > g c ; 

c ) <t 2 7^ 0, a 1 = 0. A solution of this form satisfies the equations 

1 f du G (iu) 



Qo + VlF J 2tt -iu + ef + \N nodes \a 2 \ 2 Go{iuj) 
Qf f du 1 



N c J 2tt -iu + ef + \N nodes \a 2 \ 2 Go{iu)) 
The action of this solution is given by 



(5.14) 



i s . 



-N c [ ^ log (-iu + t f + N nodes ^-G (iu) ) + Nc ^ 2 ^ - Qf e f (5.15) 
J 2n \ 2 J 9o + V9 



For a solution of this type, <Ti and a 2 are determined by a\ through the relations 



{9i - m) , 
<7\ = a 2 I = and a 2 = — a 2 



^ + 93 (5.16) 



This solution is only allowed above the line go + \f r g 1 = g c in the phase diagram. This 
is the Kondo or screening phase. The results of the previous sections apply to this 
phase. 

d ) o 2 7^ 0, d\ = 0. The SPE's also allow for solutions of this form. In principle such 
solutions are allowed above the line g — \f r g 1 = g c - The SPE's for this solution are 
exactly the same ones found in (c) with g + v^P replaced by g — \fg I - Clearly this 
solution competes with the solution of case (c) over a significant region of the phase 
diagram. In fact its allowed region is completely included within the allowed region 
for (c). However, we argue that this solution is always metastable and never occurs. 
In fact, for any finite value of y/g 2 , the interaction between the impurity and the band 
electrons will be stronger in case (c) since go + y/g 2 > go — y/g 2 , favoring the first case 
in much the same way in which the bonding state is preferred to the antibonding one 
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when the degeneracy in a two-level system is lifted by a perturbation. It is easy to 
check that this is indeed correct by looking for solutions of both types in the vicinity 
of the phase transition, for g close to g c and a/<F small. In addition, since there is 
no actual symmetry change between these two "phases", we do not expect a phase 
transition. 

We conclude that the generalized impurity model has just two phases: a phase with an 
unscreened impurity and a phase with a Kondo effect. The phase boundary is at the line 
9o + VW — 9c ■ The physics of the Kondo phase is almost identical to what was described 
in the previous sections. In fact, at the level of the N — > oo theory, the physical observables 
of the generalized impurity model can be calculated using the formulas for the single node 
model of the previous sections. The main difference is that, since the U(2) flavor symmetry 
at each node is now broken by the explicit form of the impurity Hamiltonian, operators 
with non-trivial matrix elements in that sector are now allowed. In particular, there will 
be non-vanishing, finite, off-diagonal matrix elements of the band fermion T-matrix. In 
other words, the impurity will mix band fermion states with different angular momentum. 
In contrast, processes which mix different nodes are still strictly forbidden. Finally we 
note that scattering processes of the type described here effectively reduce the number of 
independent channels. If it were not for the existence of several nodes, we would expect 
to find a Kondo effect in this phase with a completely screened impurity. Once again, the 
existence of additional channels associated with the multinode structure turns this into a 
multichannel Kondo system if inter-node scattering processes are not allowed. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have considered the problem of flux-phase fermions coupled locally to 
a single magnetic impurity. This is formally the same as the problem of fermions with a 
relativistic dispersion in two space dimensions couped to a local spin. We derived explicitely 
an effective theory in one-space dimension. 
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The physics of the effective one- dimensional theory for this problem differs from the 
conventional radial picture of the Kondo problem in several important ways. Firstly, it is 
always a multichannel system. The channels reflect the spinor structure of the nodes and 
the multiplicity of nodes. Secondly, the relativistic dispersion implies a density of states 
which vanishes linearly at the Fermi energy. Consequently, unlike the conventional Kondo 
model, the effective one-dimensional problem has a non-local coupling between the effective 
right movers and the impurity. This is the feature that drives the Kondo effect away from 
marginality and it is responsible for the phase transition between an unscreened impurity 
phase and a phase with a Kondo effect. This is consistent with earlier results of Withoff and 
Fradkin. 

However, for the case of the linearly vanishing density of states that we discussed in this 
paper, the nature of the scaling in the vicinity of this zero-temperature phase transition is 
drastically changed. We find that all physical quantities exhibit very simple scaling laws 
modified by logarithmic corrections. We found this behavior in the Kondo scale and in the 
impurity spin suceptibility. This behavior is strongly reminscent of critical phenomena at an 
upper critical dimension. One important consequence of the logarithmic corrections is that 
there are more dynamical scales and that all physical quantities are no longer controlled by 
the Kondo scale Tk alone. This is particularly clear if one compares the amplitude of the 
T-matrix and the position and width of the Kondo resonance. 

The model we studied and solved here using large- N methods is an interesting problem 
on its own right. In a subsequent publication we will report on a study of a similar model 
for magnetic impurities in <i-wave superconductors where we will draw heavily on the ideas 
that we developed here. 

It is interesting to compare the effective one-dimensional model that we derived here 
with the conventional one-dimensional models for the conventional Kondo problem. The 
standard Kondo problem is equivalent to a model in one space dimension with a single 
right mover which interacts locally with the impurity spin through the fermion spin density. 
This coupling through a density is crucial for the physics of the Kondo problem to be 
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correctly described by the model of right moving fermions. The fact that the fermions 
are chiral (namely, only right movers are present) means that , up to a Fermi velocity, 
the fermion density and current are the same observable. This model can be described 
entirely in terms of a conserved current. This is the starting point of the Conformal Field 
Theory approach of Affleck and Ludwig |21||. However, it is also crucial for the success 



of the approach of Anderson, Yuval and Hamman ||. In fact, the equivalence that exists 
between Kondo systems and problems of Macroscopic Quantum Coherence [^3| rely heavily 
on bosonization of models of fermion coupled locally to impurities through a density. The 
fact that this density is associated to a conserved current means that it cannot acquire 
anomalous dimensions. Thus, the fixed points of the conventional Kondo contains only 
marginal operators which are made marginally relevant by quantum fluctuations. 

From this analysis it is clear that it is not possible to describe the phase transition that 
we discuss here in terms of a Conformal Field Theory coupled to local boundary operators 
representing the impurity. The non-locality of the effective one-dimensional theory is es- 
sential. It is because the model is non-local that the operator that couples to the impurity 
can (and does) acquire an anomalous dimension. This is the mechanism which drives the 
phase transition. These models are not equivalent to any standard Macroscopic Quantum 
Tunneling (MQT) model. It may appear that, because the density of states vansihes like a 
power of the energy, these models could be related to a subohmic MQT system which are 
known not to have phase transitions. However, subohmic MQT models of quantum impuri- 
ties coupled to a macroscopic bosonic system, with subohmic spectral density. The models 
that we discussed here are fermionic and are not equivalent to a subohmic bosonic theory. 
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APPENDIX A: 



In section [TTT] and from Eq. ( j3.8|) we have 



i r~ 1/A v 

h = - dx — (Al) 

7T JO ^ + j2 



and 



lr 00 z/ 
J 2 = - / da; n (A2) 

In the region of integration where x » e _1 ' A , I2 can be approximated by 

J 2 » arctan (-e~ 1/A ") (A3) 



2 7T \V 

On the other hand, when we are in the regime x « e _1 / A 



v r e ~ 1/A dx 

71 A z Jo E.\ _l_^2 1™-2 



A 



+ x 2 log X 



The following inequality holds 



1 2 T-2/A < ^> A ) < - e ^ /A ( A5 ) 
7T + e 7TV 

We will be interested in the following cases. 

If v » e" 1/A , we have R(u,A) -> — e" 1/A « 1. (A6) 

If 1/ « e" 1/A , we have R(v,A) -> - e 1/A << 1. (A7) 

7T 



For 1 >> ^ >> e 1 / A , we can expand the denominator in Eq.(A"4) to get 



%A) - -Ue"^ - -^f-V e - 3/A + ...l (A8) 



7TZ/ 1 A 2 V z/ 

We can now get the behavior of Eq. (|3.1| ) in these two different regimes: for e _1 / A << v « 1, 

« e" 3/A and 1/ « - ^ » e~ 1/A if -»• -. A9 

A c 2 7tz/ 3 Vtt/ , iV 7 ^ c 2 v 7 
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the other limit corresponds to v « e 1 / A << 1; consequently the contribution from I 2 is 
neglegible and 

Of « V - e VA or ^ ^ Tr^e-VA « e -VA if % « I (A10) 
Now we want to work out similar approximations for the other SPE. We have 
= - f 1/A dx^(l-A l ° gX ^ 



1 — X 2 \ 1 — X 2 



u 2 +x 2 (i-A^y 



(All) 



in which we made use of the fact that over the integration interval, the variable x is very 
small, so that the approximations 1 — x 2 ~ 1 and — Alogx >> 1 are consistent. For the 
other portion of the integral (i.e., for I' 2 ) we are in a situation where x > e^ 1 ^ and then 
— A < 1, so we neglect the logaritmic term against 1. 
We need here the following result 

r°° X 2 log X 1 1 1 ( 7T 2 



/ dx n = n \ — + log V \ A12 

Jo 1-x 2 u 2 + x 2 2 l + i/ 2 I 2 & ) y J 



Then 



r dx — 

r°° x 2 logx 1 



I' = - dx —2— 1 — A 



logx / log a; \ x 2 



r°° , 
/ ax — 

ie-VA 1 



•• , - ! - X 2 // 2 + X 2 



l :i 

2 



T^{T + —"}-- 1/A ( 1 + i) 

-vlogv arctan (^e~ 1/A ) — i/ F (^e~ 1/A ) ( A1 3) 



The function F(z) verifies that, for \z\ > l,F(z) = F ( -M and for |z| < 1 is given by 



m = EC- 1 )" 



i 



n=0 2 ™ + 1 



2n+l 



log Z - 



2n + 1 



^ z(\ogz-l)-^z 3 (\ogz-^j+... (A14) 

Getting everything together, in the limit v ^> e _1 / A and v, A << 1, where 
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R(v, A) -> — e" 



-1/A. 



one can write the SPE in the form 

1 

go 



7T 

T 



7T 



ulogu 



which can be solved by iteration (see Eq. (|3.14j ) in section |Hl 
corresponds to v e x / A << 1, with v, A « 1. Here 



R (z/, A) 



7/ 



7T 



1/A 



The SPE now is 



(A15) 



The other limiting case 



(A16) 



(A17) 



4 £o 

since in this case, ue 1 ^ m n^f. The cubic term can be dropped rightaway. As for the other 



terms, as we approach the transition, ir^f remains fixed while A 



0. Therefore, the 



leading term is going to be the one with the factor This immediately gives the implicit 
expression for A in terms of the impurity occupancy jf- and the distance to the critical 
point ^ — ^ that we used in Eq. ( |3.9|) in section pi] . 
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FIG. 1. Phase diagram of the generalized impurity model. 
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